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O , Abstract 

I , We begin with the time-dependent electric and magnetic dipole solu- 

J^P^' tion of Maxwell's equations in Minkowski space. This Maxwell field is 

then used to determine the behavior of the gravitational field (the Weyl 
tensor) as a second-order perturbation off of the Minkowski background. 
From the Weyl tensor we go on and find the spin-coefficients and the full 
^ ' metric in this approximation. The physical meaning of many of the rela- 

te , tions is discussed. In particular we can identify the conservation law of 

^ ' ' angular momentum that contains an angular momentum fiux term. 

cn 

C^ ■ 1 Introduction 

O. 

Q^ ■ Recently Bramson[T] has investigated, perturbatively, the gravitational effects 

^^ ' due to a radiating pure electric dipole Maxwell field as source. He began 

with the pure Maxwell field in a Minkowski background space-time. He then 
k>( \ integrated the Bianchi Identities (also on the Minkowski background) for the 

j_j ■ linearized Weyl tensor with the given Maxwell field inserted into the stress tensor 

C^ I acting as the source for the Weyl tensor. He found, as expected, that indeed the 

Maxwell field did "stimulate" a gravitational response including energy loss and 
gravitational radiation. He stopped short here and made no attempt to find the 
perturbed metric. All of Bramson's work, as well as the present work, is done 
in the spin-coefficient formalism 

The purpose of the present note is three-fold: 1. we give a mild generaliza- 
tion of Bramson's work by including both electric and magnetic dipoles in the 
calculations, 2. we correct a minor error in Bramson's calculation, and 3. (our 
major contribution), we continue the perturbation calculation to include the in- 
tegration of the perturbed spin-coefficients and metric variables, thus obtaining 
the perturbed metric. 



The strategy is to integrate the Maxwell equations and the Bianchi Identities 
(as did Bramson) and then use these results in the integration of the complete set 
of spin-coefRcient and metric equations. The metric is then easily constructed. 

In Section 2 we write the time dependent electric and magnetic dipole so- 
lution of the Maxwell field, assuming no Coulomb charge. In Section 3, we 
integrate the Bianchi identities in order to obtain the behavior of the Weyl ten- 
sor (i.e., the gravitational field). In Section 4 the spin coefficients and metric 
variables are found. As will be seen, all the gravitational information will be 
coded into the Bondi shear, which is found as a function of the electromag- 
netic dipole field. In Section 5, a reality condition (an algebraic Weyl tensor 
symmetry) is applied to the known Weyl components. This leads to physical 
meaning for many of the Weyl components as well as an understanding of the 
gravitational dynamics. In Section 6 the space-time metric is easily calculated 
and given. Conclusions and comments are given in Section 7. 

2 The Maxwell Field 

In fiat Minkowski space, we choose the well-known Bondi coordinate system 
{u, r, C, (), i.e., the coordinate system where u labels the light cones with apex on 
a time-like world-line, ( = e*"^ cot(6'/2), the complex stereographic coordinates 
labeHng the individual null geodesies, while r is the affine parameter along the 
null geodesies. The Bondi null tetrad {P, n", m°, m"} is such that 1°' is tangent 
to the null geodesies. This, with a further choice of the tetrad, fixes the fiat 
space spin-coefficients (SC) as: 



(T = e = 7r = T = A = 7 = :^ = 0, (1) 



r , 2r r ^ '^ 2r r ' 



M= ^ 



In this tetrad and coordinate system, the retarded fiat-space Maxwell field 
for electric and magnetic dipole radiation can be written as[2l[3]: 



io - -^ - ^3-^ii: (2a) 



2 3 — 2 1^ 3 I'i' 



i^ = -2D'"Y^-'^. 



as 



The complex three-vector D* (the complex dipole moment) can be written 

D^=D^E+^Dli, (3) 



where D''^ and D\.j are the electric and magnetic dipole moments respectively. 
D* is an arbitrary function of the retarded time, Ur = v2w- 

For our perturbation calculations we will treat, D^ (or ^*) as first order. 

Remark: OftenfH HI [6] we find it useful to write I?' = q^^, where C is 
treated as a complex position vector. 

Remark: In the preceding as well as in what follows, we will denote the 
derivative with respect to the Bondi time as du{) = ('), while differentiation 
with respect to the retarded Bondi time is 9„^() = ('). Note that () — %/2('). 
Later in Section 5, just to view the equations with appropriate units (i.e., to 
restore the c), we make the transition (') => c^^(') in a few equations. 

3 The Gravitational Field-the Weyl Tensor 

We now turn to the linear Bianchi identities in the SC formalism with the flat 
space set of coefficients given in (1). The radial equations can be written as[3]: 
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In these equations, the fiat space operator A is given by: A = 9„ — dr and 
the constant k = 2G/c'* . If Eqs. l(2a|) had been used here a large simplification 
would have occurred. 



Now, as we are interested only in those gravitational effects which arise due 
to the Maxwell field, Eq. l|2ap . we can follow Bramson[l| in setting i/'o = 0. In 
this case, integrating Eqs.lHl), with frequent use of Clebsch-Gordon expansions[7] 
to re-express the quadratic terms, produces: 

^1 = 7r-*^(^ ^6 )^^jkY,k (5) 



24DW1' + {4 + V2)(DW^y V2DW\ „ 

+ik( ^ ^ je^jfe^ifc 

AD^'&' 30V2D'&' + (1 + 2V2){D'D^y lODW^ ^ 



v|_ evi^ ev? 
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The radially independent 'functions of integration' V'l, V'2; V's; and ■!/'4 are 
now determined by the non-radial Bianchi identities. These take the form[3]: 

5V? = 3fc0g02, (9) 

^1 = -5^-2 + 2fc</'?02, (10) 



^2 = -3^3 + kcj}24'l (11) 



■0 --/,o 



^3 = -3^4- (12) 

A straightforward 'integration' of ([9]), i.e., just the comparison of spherical 



harmonic coefficients, yields 



ij° = 3kD^&"Y'+a%\. (13) 



Here, the radially independent complex 3- vector a^{ur) emerges as another 
'function of integration', the behavior of which will be determined shortly. Con- 
tinuing on to (fTOl . we can again , by comparing spherical harmonic coefficients 
and the use of a Clebsch-Gordon expansion of the Maxwell field term, obtain: 

^2 ^ X+i^a'^' --^kiD^' D'''Ujk)Y?,, + V2k{\{D^&'')' +^-D' D^ (14) 

The scalar valued xi'^r) enters here in the same manner as the a*^ entered 
in the previous 

integration. Now, from ifTTj) . noting that since ■03 has spin-wt s = — 1 and hence 
it cannot have 
an Z = harmonic contribution, leaves us with the following relations: 

i:l = [2V2ki{D"&")'e,,k + V2ki&"Di"e,jk~a''"]Y^k (15) 

+ k[-D"'D^" - {D'D^")" - -{D"D^")']Y^^l, 
3 3 -^ 

Ah 

X = ^X'= ^&"D^"5,,. (16) 

It should be noted that the evolution equation for x is equivalent to that 
given by Bramson [I]. From a physical viewpoint it is (up to multiplicative 
constants) exactly the Bondi energy loss theorem. It coincides exactly with the 
classical E&M electric and magnetic dipole energy loss. 

Finally from Eq. lfT2|) . recalling that -04 has spin weight s — —2, it follows 
that there cannot be an / = 1 harmonic contribution in that equation. This 
allows us to determine a'' (up to initial values) as well as obtain the desired 
expression for ip1: 



04-\/2fc 



{D'D^")"' - l{D'"D^"y + -{D"D3"y 



Y2i!^ (17) 
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= V2kt{2{D"D^")" + {D'" D'")'}e,jk- (18) 



Remark: The Weyl tensor components 0]^, ■02; and -03 all appear to contain 
the I = 1 spherical harmonic contributions from the vector a^\ however 



Eq. (fT2l) puts a condition, Eq. lfTSl) . on the a'^ that ehminates the I — 1 part of 
■03. Bramson[l], improperly set the a'^ = 0, apparently overlooking this 
condition. 

Remark: There is a further equation (a reality condition) on V'2 that plays 
a very important role. It however involves a further variable (the spin-coefficient 
a'^, i.e., the Bondi shear) for its description. We postpone a discussion of this 
condition until later. 



4 The Spin Coefficients 

In order to determine the metric of the perturbed space-time with the Weyl ten- 
sor given in Section 3, we need the full set of spin coefficients. We integrate the 
spin-coefficient equations on a flat background but 'driven' by the just obtained 
Weyl tensor. Afterwards we integrate the 'metric' equations which then lead 
to the metric expression. 

Since we considered the Maxwell field to be small (first order) and since they 
entered the Bianchi Identities quadratically we must consider the Weyl tensor 
to be second order. We thus will calculate the spin-coefficients with a second 
order correction to the zeroth order flat space coefficients. 

The only SCs which explicitly remain zero in the perturbation are (k, 7r,e) 
due to our choice of the tetrad[8]. This tetrad choice also continues to fix p = p 
and T = a + (3. 

We first integrate the "optical equations" 

_=p +CTa+— ^, (19) 

^ = 2pa. (20) 

or 

for the spin coefficients p and a, the divergence and shear of /".Writing p = 
r~^ + p and cr = s, we obtain: 

p = _i+ P|_^_^^, (21) 



^=7^- (22) 

However, as the a° (our function of integration) is of second order in the 
perturbation, the third term in the expression for p can be neglected. Further- 
more, the integration constant p" can be set equal to zero by a second order 
shift of the origin of the radial coordinate r, and our final expressions become: 

p= ^- !^, a= 4. (23) 



Solving the other radial and non-radial SC equations yields the following 
results, which are in the leading terms equivalent to those of the asymptotically 
flat space-time[5: 
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>^=- + ^ + k{—-^ _^ + __)y^^, (28) 
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(30) 



2D"'D^' 8V2D"Di' + V^D^'D^ + (1 - V^){D'D^')' ^ 

— ki{ T — eiitri, 

70V2D''D^ + {3-40V2)iD'D^y 11D'D\ ^^ 

The radially independent constants entering into the formulation are gov- 
erned by the following relations which emerge directly from the non-radial SC 
equations: 

^° = V2ga°', 

^0 = -V2A°' = -2a"". 

First, however, by comparing the ■04 from Eqs. ifT7|) and l|3ip . i.e.. 



V;° = -V2fc[^(5"'i?^")' - {D'&T' - \{D''&")"]Y^f, 
with i/;" = —A = — o- , we obtain the asymptotic Bondi shear 



(32) 



?" = ^^^ j{&"Dndur - (i?^5^")' - |(i?^'5-'")]r2^, (33) 

and the Bondi news function, 

^' = ^mJ(^'"^^") - (i?^^^")" - ^(i^^'^^")']r2^. (34) 

These are basic and important relations. 

With these results all the functions of integration are now known in terms 
of the Maxwell field, or equivalently in terms of the complex D^ and its time 
derivatives. 



5 The Reality Condition 

In the previous section we mentioned that there existed a further condition on 
the asymptotic Weyl tensor, namely a reality condition. By defining the mass 
aspect, ^, by 



* = V-" + 3V + V2a"{a°y, (35) 

the reality condition is simply 

* = *. (36) 

Before we find the consequences of this condition it is worthwhile to see the 
physical meaning of the mass aspect. 

When '^ is expanded in spherical harmonics, 

q> = ^i>° + ^S'Y^^ + -i''^Y^^^+..., (37) 

one has, from the work of Hermann Bondi[9l [3j, the physical meaning of the 
I = (0, 1) harmonics as the mass and the linear 3-momentum: i.e., 

*° = -^4^M, (38) 

^, ^ -^P\ (39) 

The first result from the reality is that the M and P' are real. We will see 
shortly that it also leads to an expression for conservation of angular momentum 
in the form of a flux law. 

Returning to the expression for ^, we see that since a'^ is 2"'' order the 
termV2cr''(CT")' can be ignored and the condition becomes 

Substituting the expressions for ip2 and cr° (i.e., Eqs. (fT4| and l|33p ) into ^ we 
obtain, after some cancellation of terms 

* = X + (^a" - 2kzD^'&"c,,k)Y,", + k^ j {&" UndurY^,^ . (40) 

The I = 2 harmonic is real without any added conditions. This follows from 
the symmetry of 1^2*^^ • 

The I = harmonic condition is just 

*°-X = X = -^M, (41) 

while the I = 1 condition, with a^ = a^ + ia'}, is that 

af = 2V2k[D%D%' + DtjDi'i]e^,k ■ (42) 

Comment: In an earlier section we saw that the complex a'^ was de- 
termined by £^(j>. . (fTSjl ■ i.e., by 

a'="' = 2V2ki{D''D^")"e,jk + V2ki{D'" D^")'e,jk. (43) 



It was a relief to us to see that the imaginary part of Eq. ^AS^ coincided 
with the reality condition, Eq.^^. For completeness we give the real part of 
Eq.^: 



,kiir 



= 2V2k{D'l,D%y"e,,k - V2k{D%'D^;;ye 



Ml ^ijk- 



The linear momentum given by 



*' = -^P"" = [^a% - 2k{D^;,D%)e,,,]'. 



(44) 



(45) 



The last piece of information that we obtain from the reality condition comes 
from the conventional identification in linear theory of source angular momen- 
tum as observed from infinity, namely: 



J'^ 



12G 



■V-ifcl/- 



(46) 



From Eq. (|T3|) we see that 



fc _ V2^ , 



r ^ 



12G 



*7: 



(47) 



so that, from Eq. l(42|) . with the restored c's from (') => c ^('), we obtain 



2 
3c' 



J'=' = :^[i)j^< + D^<']e,,fe. 



(48) 



In other words the reality condition appHed to the I = 1 harmonic, i.e., 
Eq. (|42l) . describes the loss of angular momentum. This agrees exactly with the 
classical loss of angular momentum due to electromagnetic dipole radiation\10\. 

6 The Metric 

In our Bondi coordinate system, the inverse metric takes the formjS]: 
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where all the metric components are given in terms of metric coefficients U, uj, 
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2iU-uju), 



^^^X^-iQ^^+iut), 



^AB 



-(e^f +^^C-'), 



(49) 



These metric coefficients satisfy a set of radial and non-radial equations given 
in terms of the spin coefficients just calculated O [8]: 
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-g:^ - puj + auj - {a + (3), 



(50) 



<SX^-A^^-(M + 7-7)r + Ar, 

<5e^-je^-(^-«)e^ + (a-/3)^^, (51) 

du — Suj ~ {P — a)u! + (a — /3)w + fi — fi, 
SU — Auj = {p + J — j)lli + Xlu — D. 

Still using the flat-space derivative operators, these equations can be inte- 
grated as was done earlier, we obtain (using certain standard coordinate condi- 
tions for the sphere metric) the metric coefficients: 

^oc = 0, f '^ = -P, (52) 

^^(1 + CC), 

^ = ^^ ;^' (53) 



w 



So-" ^Z-? , ADW3' V2DW\ , , , 



^■^ ^ r[|i + ..(^-^lr)^*n^' (-) 



Lg^s ^ 3J.4 lOrS 

-K — ^-3 T7r:r)^2, 



C/=-l-^^+A (56) 

2r 

where we have bundled higher powers of r into the term A given explicitly by: 
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^ = —^ o(^ —?, + 1.^4 )^n (57) 



6r^ 3 Sr^ 

ki{2D3'D' - (DWy) 
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We can now compute the inverse metric to second-order in the perturbation, 
where the components given in Eq. l|49p take the form: 



(58) 



," = -2-M±i^ + 2A 
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Calculating the inverse of g .i.e 



gab 



.900 1 goc 5oc 
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.9oc 5^^ g^-j 

5oc 5cc ffcc 
we obtain the metric, now fully expressed in terms of Maxwell field or in terms 

of the complex dipole, D* : 



+r^ik[—{D'D^" + D'D^ 



V2{D^'Dn'^o 



500 = 2 + r-l{2x+[^(fl" + a'=')-2fc^((i?^'^^■')')e.J■fc]n°J (62) 



]>?«, - 2 A 
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,C0 = -P-'i^'r^'^'-^Y^^^ki^^'-y^)e^,,Y^ (63) 



1 

3^2 ISr'^ 



+fc( — : ^3 — + ^7i?^)^4), 
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2DW^" 4:V2DW^' ADW^ . 
3j,2 15r'^ 



-^-(5*"+ 4^V,' + .«i^-i^;ii^)e«.y- (64) 

+ fc(^ 512 + ^71^)^27, ) 
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We recall that 



a"= ^[^j D^"D^"du,-{D^Dn'' li&'Dn]Y,%. (66) 



Eqs. l(62l) - l(66|) describe the metric resulting from the perturbation off flat 
space induced by the electric and magnetic dipoles Maxwell fleld. 

7 Discussion and Conclusion 

In the present work we have extended ideas and work of Bramson. Starting with 
the radiating dipole Maxwell field, we found in addition to the 2nd order Weyl 
tensor, the spin coefficients and metric up to second order. We also found the 
physical content of the I = 1 spherical harmonic of the reality condition to be 
the angular momentum conservation law. Thus, in our approximation we can 
fully characterize, up to second order, the changes imposed on the Minkowski 
space-time by the electromagnetic perturbation. 

Bramson, by making the error of setting a'^ = 0, missed the interesting 
physics contained in a'^'. Furthermore, by not continuing on to the spin-coefficient 
equations he had no means to investigate the reality conditions and would have 
missed the conservation of angular momentum. One of the most important of 
the spin-coefficients is the Bondi asymptotic shear. We found that the shear (or 
its time derivative, the Bondi news-function) is fully determined by the electric 
and magnetic dipole moments. We also saw the Bondi mass loss equation yields 
the classical energy loss from electrodynamics as well as finding the Bondi linear 
momentum also determined by the Maxwell field. 



13 



An interesting observation is that in the Bondi mass loss formula, the energy 
is lost only through the Maxwell dipole radiation and there is no term involving 
the resultant gravitational radiation, via the Bondi news-function. The reason 
for this is that the gravitational radiation is 4th order and thus not visible at our 
approximation. It however does suggest that whenever there is electromagnetic 
radiation it will always be accompanied by weak gravitational radiation. 

This approximation scheme can be applied to any number of perturbation 
calculations. In a future paper, we will apply this method to a more general 
perturbation off the Schwarzschild solution. One should be able in principle to 
do similar things from a Reissner-Nordstrom or Kerr background as well, and 
in future work we will investigate these scenarios. 
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